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Abstract 

This paper deals with fractional-order controlled systems and fractional-order controllers in the dis- 
crete domain. The mathematical description by the fractional difference equations and properties of 
these systems are presented. A practical example for modelling the fractional-order control loop is 
shown and obtained results are discussed in conclusion. 
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1. INTRODUCTION 

Fractional calculus was used for modelling of physical systems, but we can find only few works dealing 
with the application of this mathematical tool in control theory (e.g. |l], |^, |j| ||, [?], ^). These works 
used continuous mathematical models of fractional order. The fractional-order systems have a unlimited 
memory, being integer-order systems cases in which the memory is limited. 

It is necessary to realise that the case of discrete fractional-order systems is very important for 
their description to have a finite difference equation. Such an equation can be obtained by numerical 
approximation and by the inverse Z-transform of a discrete transfer function. 

The aim of this paper is to show, how by using the fractional calculus, we can obtain a more general 
fractional-order models of the controlled objects and general structure for the classical PID controller in 
the discrete domain. 

2. DISCRETE FRACTIONAL CALCULUS 

The idea of fractional calculus has been known since the development of the regular calculus, with the 
first reference probably being associated with Leibniz and L'Hospital in 1695. 

Fractional calculus is a generalisation of integration and differentiation to non-integer order funda- 
mental operator a Df, where a and t are the limits of the operation. The continuous integro-differential 
operator is defined as 

f & 5R(a) > 0, 

a D?={ 1 »(a)=0, 
{ J*(dT)- a R(a) < 0. 
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The two definitions used for the general fractional differintegral are the Griinwald-Letnikov (GL) definition 
and the Riemann-Liouville (RL) definition JgJ . The GL is given here 



a D?f(t) = ]imh- a J2(- 1 ) i 

h—>0 — ' 
3=0 



f(t-jh), 



where [x] means the integer part of x. The RL definition is given as 



aD?f(t) 



1 



d" 



f(r) 



-dr, 



(1) 



(2) 



T(n - a) dt n J a (t - T ) a - n+1 

for (n — 1 < a < n) and where T(.) is the well known Euler's Gamma function. 

In general, the approximation for fractional operator of order a can be expressed by generating 
function where is the shift operator. This generating function and its expansion determine 

both the form of the approximation and the coefficients Q . 

For discrete approximation of the time derivative, we can use the generating function corresponding 
to the Z-transform of backward difference rule, 



w(*- X ) = 



1 



T 



(3) 



and performing the power series expansion (PSE) of (I — z~ x ) ±a , we obtain the Z version of the GL 
formula by using the short memory principle |l , for the discrete equivalent of the fractional-order integro- 
diffcrential operator w(2;~ 1 ) ±Q , 



(wC*- 1 )) 



l\\±a 



T^^(-l) 

3=0 



±a 



At}-3 



where T is the sample period, L is the memory length, [x] is the integer part of x and (—1) 
a binomial coefficients cf \ (j = 0, 1, . ..). For its calculation we can use the following expression: 



f±a\ 



» 



» 



1 + (±a) \ ( Q ) 
J J °'- V 



(4) 
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(5) 



Another possibility for the discrete approximation is the use of the trapezoidal (Tustin) rule as a gen- 
erating function for PSE. For the discrete equivalent of the fractional-order integro-differential operator, 
we can write a general formula Q 

for obtaining the coefficients and the form of the approximation. 

A detailed review of the approximation methods (Calson's, Chareff's, CFE, Matsuda's, 
continuous and discrete fractional-order models was done in work fl(i|| . 



,) for 



3. FRACTIONAL-ORDER CONTROL CIRCUIT 



We will be studying the control system shown in Fig.l, where G c (z) is the controller transfer function, 
G s (z) is the controlled system transfer function, W(z) is an input, E(z) is an error, U(z) is the output 
from controller and Y(z) is the output from system 



W(z) + ^ E(z) 


G c {z) 


U(z) 


G s (z) 


Y(z) 























Figure 1: Feed - back control loop 
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3.1 Fractional-order controlled system 

The fractional-order controlled system will be represented with a fractional model with the fractional 
differential equation given by the following expression (oD^ = D^): 

a n D^y(t) + ... + ai D^y(t) + a Q D f ^y{t) = b n D^u(t) + . . . + 61 D^u(t) + b Q D?°u(t), (7) 

where ft, (k = 0, 1, 2, . . .) are generally real numbers, ft > . . . > ft > ft, a m > . . . > a% > ao and 
flfe) bk (k = 0, 1, . . .) are arbitrary constants. 

For obtaining a discrete model of the fractional-order system (^), we have to use discrete approxima- 
tions of the fractional-order integro-diffcrcntial operators and then we obtain a general expression for the 
discrete transfer function of the controlled system 

s[ ' a„(w(z- 1 ))' 3 " ■ + ... + a 1 (u(z- 1 ))Pi+a (u J (z-i))P° ' [ ) 

where (^(z^ 1 )) denotes the discrete operator, expressed as a function of the complex variable z or the 
shift operator z _1 . 

For discrete time step k, according to relation (Q) and the inverse Z-transform of the difference equation 
(Q), we can write the fractional-order difference equation in the form 



i=o 3=0 3=0 

6 m T- Q ™ £ cj *" Vi + . . . + h T- ai J2 + ^ T- ao c? 0) ^-i- (9) 

3=0 j=0 j=0 

From relation (j^) the following form of the difference equation can be obtained 

^ = ^ v*. W) ' (10) 

for k — 2, 3, . . ., where yo = and yi = 0. 

3.2 Fractional order controller 

The fractional PI X D S controller will be represented by discrete transfer function given in the following 
expression: 



G c (z) = = K + T \ + ^^(z- 1 )) 5 , ( I I ) 

E{z) (uj(z 1 )) A 



where A and S are arbitrary real numbers (A, 5 > 0), K is the proportional constant, Tj is the integration 
constant and Td is the derivative constant. 

Taking A = f and 6 = 1, we obtain a classical PPD controller. If A = and/or Tj = 0, we obtain 
a PD* 5 controller, etc. All these types of controllers are particular cases of the PI X D S controller, which is 
more flexible and gives an opportunity to better adjust the dynamical properties of the fractional-order 
control system. 

Similar to methods for controlled system, they can be used for obtaining the difference equation of 
the fractional-order controller (pT|). 



4. ILLUSTRATIVE EXAMPLE 

We give in this section an example of modelling the stable dynamical system by using fractional calculus 
in discrete domain. The fractional-order control system consists of the real controlled system with the 
coefficients: 

a 2 = 0.8, a x = 0.5, a = 1.0, ft = 2.2, ft = 0.9, ft = 0, b = 1.0, a = (12) 
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and the fractional-order PD controller, designed on the stability measure St = 2.0 and damping measure 
£ = 0.4, with the coefficients: 

K = 50.0, T d = 5.326, 6 = 1.286. (13) 
The fractional-order differential equation of the closed control loop has the form (o-D^ 1 = D 

a 2 D? 2 y{t) + ai D^y(t) + T d D s t y(t) + (a + K)y{t) = K w{t) + T d D s t w{t). (14) 

The resulting differential equation ( |l4| ) can be rewritten to the fractional difference equation of the closed 
control loop in general form via inverse Z-transform of the relation (|]): 

Kw k + T d T- s jh cfw^ - a 2 T-fc £ cf\ k _ - a x T^ £ cf^Vk-j ~ T d T~ s £ cfy^ 
_ j=o J=l £=1 i=i 

W ~ a 2 T-^ + a.T-^ + TtT-stf + (a + K) ' 

(15) 

for [k = 1,2, . . .), where yo — 0, wq — 0, = and u>fc = 1, for (k — 2, 3, . . .). The binomial coefficients 
were calculated according to relation (0) . 



5. CONCLUSION 

The above methods make it possible to model and analyse (simulate) fractional-order control systems 
in the discrete domain and also to realise digital fractional-order controllers. We have shown the sur- 
vey of approximation methods and Z-transform method as good approximations of the fractional-order 
operator a D". These methods for discretisation of fractional calculus are suitable for realisation and 
implementation of fractional-order controllers, because this controllers are more robust than classical 
one. 

Stability investigation of discrete fractional-order control system and practical simulation will be 
studied in further work. 
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